Environment-induced two-mode entanglement in quantum Brownian motion 
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The time evolution of quantum correlations of entangled two-mode continuous variable states 
is examined in single-reservoir as well as two-reservoir models, representing noisy correlated or 
uncorrelated non-Markovian quantum channels. For this purpose the model of quantum Brownian 
motion is extended. Various separability criteria for Gaussian continuous variable systems are 
applied. In both types of reservoir models moderate non-Markovian effects prolong the separability 
time scales. However, in these models the properties of the stationary state may differ. In the two- 
reservoir model the initial entanglement is completely lost and both modes are finally uncorrelated. 
In a common reservoir both modes interact indirectly via the coupling to the same bath variables. 
Below a critical bath temperature entanglement between the two modes is preserved even in the 
steady state. A separability criterion is derived, which depends on the bath temperature and the 
response function of the open quantum system. Thus, the extended quantum Brownian motion 
model of a two-mode continuous variable system in a common reservoir provides an example of 
environment-induced entanglement. 



I. INTRODUCTION 

During the last decade quantum information and 
computation has been extended from discrete systems 
to quantum systems with continuous variables such 
as position and momentum or the amplitudes of elec- 
tromagnetic field modes. This quantum information 
theory of continuous variable systems has received 
much attention in the past few years [3, 0, 0| and has 
found various applications in quantum cryptography 
and quantum teleportation 0, Q. In this context, 
Gaussian states play a prominent role since they can 
be easily created and controlled experimentally and are 
less affected by decoherence. Great advances have been 
made in characterizing the entanglement properties of 
two-mode Gaussian states by determining the necessary 
and sufficient criteria for their separability 0, 0] and 
by developing quantitative entanglement measures 
[1, Q . The prototype of these states are the two- mode 
squeezed states, which have been successfully produced 
via nonlinear parametric down conversation and applied 
to quantum teleportation [lol . [ll| . 



Due to the unavoidable interaction with the environ- 
ment, any pure quantum state used in some quantum 
information process evolves into a mixed state. Thus, 
a realistic analysis of continuous variable quantum 
channels must take decoherence and dissipation into 
account. Within the theory of open quantum systems 
[iH fl^ the dissipative dynamics are mainly described by 
master equations of the reduced density matrix. Initial 
quantum superpositions are destroyed and quantum 
correlations are lost during characteristic decoherence 
and separability time scales. The Markovian time evo- 
lution of quantum correlations of entangled two-mode 
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continuous variable states has been examined in single- 
reservoir [3, [lj| and two-reservoir models @, [H, H3) 13 > 
representing noisy correlated or uncorrelated Markovian 
quantum channels. Quantum correlations are found to 
be better preserved in a common reservoir. Additionally 
the coupling to the same bath variables might generate 
new quantum correlations between the parts of the 
subsystem. This effect of environment-induced entan- 
glement has already been studied for discrete systems 
111 [IIS, and can lead to asymptotically entangled 
states [2lT|22|, [H Hi]. The underlying Born-Markov 
Approximation assumes weak coupling between the 
system and the environment to justify a perturbative 
treatment and neglects short-time correlations between 
the system and the reservoir. This approach has been 
widely and successfully employed in the field of quantum 
optics [25[ where the characteristic time scales of the 
environmental correlations is much shorter compared to 
the internal system dynamics. Challenged by new exper- 
imental evidence a growing interest in non-Markovian 
descriptions can be observed. Very recently some 
phenomenological [26|, [2?| and microscopic models 
[HI, US US, El of non-Markovian quantum channels have 
been proposed. Using the analogy between the Hilbert 
space of quantized electromagnetic fields and the Hilbert 
space of quantum harmonic oscillators, the Caldeira 
Leggett model of quantum Brownian motion [32 . HH, [34| 
can be extended to describe the entanglement dynamics 
of two-mode squeezed states. 



In this paper, the time evolution of quantum correla- 
tions of initially entangled two-mode continuous variable 
states is numerically examined in a common reservoir 
model. The focus is on non-Markovian influences and 
strong coupling effects. The non-Markovian dynamics 
are described by an extended two-mode version of 
Hu-Paz-Zhang master equation of quantum Brownian 
motion. In the single or common reservoir model both 
oscillators (modes) are coupled to the same reservoir 
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variables, whereas in the two-reservoir model each 
oscillator (mode) is interacting with its own independent 
reservoir. In both types of reservoir models moderate 
non-Markovian effects prolong the separability time 
scales which also depend on the interaction strength 
between the system and the environment. However, 
in these models the properties of the stationary state 
may differ. In the two-reservoir model the initial 
entanglement is completely lost and both modes are 
finally uncorrelated. In a common reservoir both modes 
interact indirectly via the coupling to the same bath 
variables. Therefore, new quantum correlations may 
emerge between the two modes. Below a critical bath 
temperature entanglement is preserved even in the 
steady state. A separability criterion is derived, which 
depends on the bath temperature and the response 
function of the open quantum system. Thus, the ex- 
tended quantum Brownian motion model of a two-mode 
continuous variable system in a common reservoir 
provides an example of environment-induced quantum 
two-mode entanglement. 

The paper is organized as follows. In section II, we 
briefly review various separability criteria for two-mode 
Gaussian states and the Markovian separability times 
for two-mode squeezed states in single- and two-reservoir 
models. In section III we shortly describe the Hu-Paz- 
Zhang master equation of quantum Brownian motion 
which is the basis for studying non-Markovian effects. We 
resume the extended, two-mode-version of the Caldeira 
Leggett model for single and two-reservoir models and 
introduce a modified common-reservoir model. In sec- 
tion IV we present and discuss the numerical results of 
the entanglement dynamics of two mode squeezed states 
in the modified common-reservoir model. Different sce- 
narios are analyzed, including the case of noise-induced 
steady state entanglement. A simple separability criteria 
for the stationary two-mode state is derived. Finally, a 
brief summary is given in section V. 



II. CONTINUOUS VARIABLE SYSTEMS 

A. Separability criteria for two mode Gaussian 
states 

In the following we review the separability criteria for 
a special class of continuous variable systems - the two 
mode Gaussian states. A Gaussian two mode state with 
coordinates q\ , (72 und momenta p\ , P2 has a Gaussian 
Wigner function in semi-classical phase space 
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with X = (<Zi,pi, 92, P2) and commutation relation 
[liiPj] = ifi^ij for i,j = 1,2. It is completely charac- 
terized by its first and second moments. 



V = 



AC 



(2) 



(where all first moments (Xj) have been set to zero) con- 
tains four local symplectic invariants in form of the de- 
terminants of the block-matrices A, B, C and V. A 
Gaussian continuous variable state is separable if and 
only if the partial transpose p of its density matrix p is 
non-negative(PPT-criterion) . Based on the above invari- 
ants Simon [6j has derived a PPT-criterion for bipartite 
Gaussian continuous variable states, 

S(t) = detV-i(detA + detB + 2|detC|)+— > (3) 
4 16 

which is also a necessary separability criterion for non- 
Gaussian states. The PPT-criterion has a geometrical 
interpretation as mirror reflection of the Wigner function 
in phase space. In case of a Gaussian two mode state the 
partial transpose coincides with a change of the signs in 
those elements of the covariance matrix, which connect 
the momentum of the first mode to the coordinate of the 
second mode. Equivalent to that is the criterion > i 
for separability where v± are the symplectic eigenvalues 
of the partial transposed density matrix p: 



1 

V2 



Ay ± \/Ap - 4detV" 



(4) 



with Ay = det A + det B - 2 det C. With the symplectic 
eigenvalue V- the logarithmic negativity as a quantitative 
measure of entanglement can be defined by [H, [35| : 



Ej^(p) = max{0, — ln2£_} . 



(5) 



In order to apply the PPT-criterion the complete knowl- 
edge of all second moments is required. For practical pur- 
poses there are also weaker separability criteria based on 
linear or quadratic combinations of only a few elements 
of the covariance matrix. Duan et. al. Q derived such 
a criterion. Starting point is the definition of a pair of 
EPR-like operators u = \a\q\ + ^q-i an d v = \a\pi — ^p2 
with a £ K \ {0} fulfilling the commutation relations 
[<jkiPj] — i&kj, i,k — 1,2. Then, for every separable bi- 
partite quantum state p, the sum of the variances fulfills 
the relation 
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) P > a 2 + — 



(6) 



For Gaussian states this set of inequalities (for all possible 
values of a) completely characterizes the set of separable 
states and is equivalent to eq. ([3]). The product version 
of this inequality for a — 1 is given by [la, H3, [38| 



((Auf) p ■ ((Avf) p > 1, 



(7) 



which is a special case of the generalized curved quadratic 
entanglement witnesses (39j . In general, the product wit- 
nesses are stronger tests than the respective linear tests. 
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Thus, the set of entangled covariance matrices detected 
by this quadratic test is strictly larger than that detected 
by linear combinations of second moments [3!| . 
We also will apply a classification scheme of quantum 
states based on marginal and global purities of a bipartite 
Gaussian quantum system whereas p = Tr[p 2 ] is the pu- 



rity of the total system and p\ = Tr[p 2 ], p2 = ir^J are 
the purities of the reduced density matrices p\ — T12 [p] , 
Pi = Tri[p] of the two- mode system. According to ref. 
[40 , l4l| , a two- mode Gaussian mixed state is separable if 
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and it is entangled for the case that 
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< P < 



P1P2 



P1P2 - I Mi - M2I 



In between there is a coexistence region 
M1M2 „ . M1M2 



Ml + M2 - M1M2 



< p < 



m! 
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(8) 



(9) 



(10) 



where both separable and entangled states can be found. 
The calculation of the global and marginal purities thus 
provides analytical bounds on the entanglement of Gaus- 
sian states. 



B. Markovian separability times for two-mode 
squeezed states 



with squeezing parameter £, the initial covariance matrix 
is already given in standard form with elements a = b = 
\ cosh(2|£|) and c± = ±i sinh(2|£|). Applying the above 
separability criterion, a Markovian separability time can 
be derived from the separability function ([3]) which is 
reduced to S(t) = e^^e" 7 * + A(t) - 1. The time n 
after which the initial entanglement is lost is given by 
the condition S(ti) = and reads 



Ti = — In 

7 
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1 - e ~ 2 l*l 
2n 



(14) 



0, EE E3 for the two-reservoir model yielding ri — > 00 
for n — > 0. In the case of a single reservoir model where 
the corresponding master equation contains additional 
terms, the time evolution of the covariance matrix ele- 
ments is slightly different and results in the separability 
time [ill 



T2 



2 7 V 2n 
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1 



»2|f| 



(15) 



Thus a common reservoir extends the Markovian sepa- 
rability time. Furthermore the initial entanglement is 
partially preserved even for n > if the initial squeezing 
exceeds a critical value £ c = | In (2n +1). These results 
change if the dynamics are described without the Born- 
Markov assumption. 



The separability criterion (|3|) is simplified if the covari- 
ance matrix V is transformed to standardform by sym- 
plcctic transformations: 



fa c + 

a c_ 

c + 6 

V c_ 6 



(11) 



where the elements a, 6, c+ and c_ of V^t are determined 
by the four symplectic invariants det A = a det B = b, 
detV = (ab — c+)(ab — c?_) and det C = c + c_. In a 
Markovian two-reservoir model the dynamics of the two- 
mode state are described by the quantum optical master 
equation of the damped quantum oscillator (interaction 
picture) 



^ \ (nL[a]]p +(n + l)L[ aj ]p 



3=1 



(12) 



III. NON-MARKOVIAN ENTANGLEMENT 
DYNAMICS 

A. HPZ master equation of quantum Brownian 
motion 



Non-Markovian effects are discussed here on the ba- 
sis of the Caldeira-Leggett model of quantum Brown- 
ian motion [32|, l42Ll43Tl o ften referred to as independent- 
oscillator-model [33[~ 44 1 . It is a system plus reservoir 
model where the total Hamiltonian consists of three parts 



H — H s + Hb + H 



int • 



(16) 



with H s as Hamiltonian of the subsystem which interacts 
via the Hamiltonian H lnt with a bath that is described 
by a collection of a large number of harmonic oscillators 
Hi, = foojifi^b + 1). In detail the Hamiltonian of the 
Caldeira Leggett model is given by 



with mean bosonic occupation number n = (e@ huJo — 
damping constant 7 and I/[o]p = lopcfi — <yop — po^o. 
The time evolution of the matrix elements is then given 
by a(t) = ae-T* + A(t), b(t) = be" 7 ' + A(t) and c±(t) = 
±c±e-* )t with A(i) = (2n+ 1)(1 - e" 7 *) In the case 
of a two mode squeezed vacuum state 



|0)i|0) s 



(13) 



2 N 



Pi , 

2m, 2 



Ciq 



(17) 



where q and p are the Heisenberg-operators of coordi- 
nate and momenta of the Brownian oscillator moving in 
an harmonic potential V{q) 



hmu>Qq 2 and coupled to 
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a bath of N independent harmonic oscillators with vari- 
ables Xj, Pi and frequencies u>i. The bath is characterized 
by its spectral density 



= 7T 



JV 9 

2muj. 

i=i 



<"-"i) = ^?F2- ( 18 ) 



The interaction is bilinear in the coordinates q and Xi of 
the subsystem and the bath respectively with coupling 
parameters Cj. The self-interaction term (proportional 
to q 2 ) in the Hamiltonian 



E 



2m l ujf 



(19) 



renormalizes the oscillator potential to ensure that the 
observable frequency is close to bare oscillator frequency 
luq. From influence functional path integral techniques 
Hu, Paz, Zhang have derived the master equation 



P = 



2ih 



[1 2 ,p] 



2ih 



.1, {p, p}] 



qp{t) [q,\p, P ]]- " 



n 2 



Dp(t) 

h 2 



k, [q,p]], 



(20) 



with [ . ] and { . } denoting commutator and anti- 
commutator respectively. This master equation is valid 
for arbitrary coupling and temperature. The non- 
Markovian character is contained in the time-dependent 
coefficients which read in expansion up to the second or- 
der in the system-bath coupling constant (34] |: 

7p(*) = IT— I dt'L(t') sinuot', (21) 

Sn 2 (t) = — - / dt'Lft') cos oj t', (22) 
m rvm J 



D qP {t) 



moj j 
ft 



dt'K(t')smu t' : 



D p (t) = ( dt' K{t') cos uj t' , 
Jo 



(23) 
(24) 



where L(t) = i([r,(t), T)(0)]) and K(t) = £<{rj(t),»K0)}) 
are connected to the spectral density (fT5)) by 



h f°° 

L(t) = — I duj J(u>) sinwt, 



(25) 



o 



h f°° l 

K(t) = - / duo J(ui)coth(-f3huj)cosujt. (26) 
n Jo 2 

K{t) is the correlation function of the quantum noise 
term r] resulting from averaging over the initial thermal 
bath distribution. The exact expressions of the HPZ- 
coefficients are related to the Green's functions of the 
corresponding quantum Langevin equations (45l . |46| . The 
entanglement properties of the joint state of the oscillator 
and its environment have been studied in ref. 14711 . 



B. Two- reservoir model 

The dynamics of two identical, not directly interacting 
modes (with coordinates and momenta qj,Pj, j = 1, 2) in 
two uncorrelated reservoirs is modeled by the interaction 
Hamiltonian 



-<?i 



OO 

E 



<?2 E 



(27) 



with (x^x^ + xfx b 



= 



, .. , Vi,j. The master equation of 

the reduced density matrix is then given by the sum of 
the master equations of two single modes [28| : 




Qj>{Pj>P}] 



(28) 



The time dependent coefficients are given by 7 ? (t) = 
uj 2 + Sfl 2 (t) - 7r/m and eq. ([21]) to (|24"|). The time evo- 
lution of a two-mode squeezed state in two uncorrelated 
non-Markovian channels has been studied recently in ref. 
[28| . The authors derived a non-Markovian separability 
function which shows oscillations in case of an artificial 
out of resonance bath with T -C loo- In this two- reservoir 
model the initial entanglement is always completely lost 
and both modes are finally uncorrelated (even at zero 
temperature while t\ — > oo). This may not be the case 
in a common reservoir model as will be shown in the next 
sections. 



C. Common-reservoir model 

The dynamics of two identical, not directly interacting 
modes in a common reservoir is modeled by the interac- 
tion Hamiltonian 



Hint = -glE 



Ciqi 



oc 
i=l 



(29) 



The corresponding master equation for a two-mode sys- 
tem with H s = H s i + H S 2 is then given by [48| 



P = -r [H s ,p] 
in 



ft 2 



M5Q 2 (t) 
ih 

[R, [Pr,p]}- 



[R 2 ,p} + 
h 2 



[R,{Pr,p}} 



ift 

{R,{RiP}} (30) 



with 5Vl 2 (t) = 5fl 2 (t) — -fT/m. As one can see immedi- 
ately, the master equation (|30|) contains only the normal 
coordinates R = (q± + q<i)l1 and Pr = pi + P2 but not 
the relativ coordinate x — qi — q2- This results from the 
specific coupling ~ Xi{q\ +52) and from the choice of two 
identical modes with uj\ — uj2- In this case, the motion of 
the relativ coordinate is not accompanied by dissipation 
[49l . [501 ] . Additional assumptions are therefore necessary 
to describe the relaxation dynamics completely. 
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D. Modified common-reservoir model 



The Hamiltonian (|29|) describes a bilinear coupling in 
the coordinates of the two modes and the coordinates of 
the bath variables. The subsystem is thus effectively cou- 
pled by the center-of-mass R = 5(91+92)- The motion of 
the relative coordinate however is unitary. Therefore, in 
the case of an initial two mode squeezed state only half 
of the entanglement is lost in the stationary state regard- 
less of the coupling strength [3l[. In order to be able to 
model a complete loss of the initial entanglement, an ad- 
ditional assumption about the dynamics of the relative 
coordinate is necessary if the two modes are identical. To 
be able to model the dissipative dynamics for arbitrary 
coupling we choose an interaction Hamiltonian similar 
to the original Caldeira-Leggett model substituting the 
variable q by R. A renormalization term proportional to 
R 2 ensures that the frequency shift is compensated in the 
stationary state. The modified interaction Hamiltonian 
thus reads 



where x(uj) — \Mu> 2 — Moj — iuij(u>)] 1 with 7(0;) = 
J °° dt"i{t)e %ult is the susceptibility of the non-Markovian 
damped harmonic oscillator. The motion of the relative 
coordinate x = (91—92) is decoupled from the motion of 
R. In order to study the case of two identical modes with 
dissipation in the relative motion additional assumptions 
make sense. It seems plausible to assume a weakly dissi- 
pative dynamics of the relative coordinate given by the 
solution of the Born-Markovian master equation of the 
damped harmonic oscillator in the form 

( x 2 (t)) = (xDe-^ + 2(1 - e-^*)(2n + 1), (37) 



where 7 P = lim t 



7 P (i) is the stationary (Markovian) 



value of the coefficient (1211) . The initial value 



) m case 

of a squeezed two mode state is given by (xq) — e -2 ^. 
The dynamics of the normal coordinates is uncorrelated 
with ({x,R}(t)) — 0. Of course, correlations in the orig- 
inal coordinates are still present. 



2J CiXiR 



00 r 2 



:R Z 



(31) 



IV. NUMERICAL ANALYSIS OF THE 
MODIFIED COMMON-RESERVOIR MODEL 



to guarantee that the eigenvalue equation for the eigen- 
frequencies v of the total system 

[v 2 -ujf + h{v)] [v 2 -lo 2 + h(u)] - h 2 {v) = (32) 



with h{v) = J2i 



has just eigenvalues v > 



(all masses m, have been set to one). The frequency 
shift in the stationary state is compensated and the ob- 
servable and the bare oscillator frequency coincide. The 
Hamiltonian of the total system consists of the Hamil- 
tonian of the two modes H s = H s i + H s2 , the Hamil- 
tonian Hi, = J^i faiJi(b>b + 1) of the bath modes as well 
as the interaction Hamiltonian (|3Tj) . The system of cou- 
pled Heisenberg equations of motion can be solved equiv- 
alently to the case of a single Brownian oscillator and 
leads to two coupled quantum Langevin equations for 
the two modes j = 1,2: 



Qj = -UjQj U / dt 77 + — > ( 33 ) 



M 



1 , 7 (t- t')R{t') , 

M 



M 



with 7(4) = 7re _rt . In the case of two identical modes 
with u>j = ojq for j — 1,2 and total mass M = 2m the 
system of coupled Langevin equations (|33|) reduces to a 
quantum Langevin equation of the center-of-mass motion 

R -~"° R Tl X dt M + 1*' (34) 



with stationary correlations 
(R 2 ) 



h f°° 1 

— / duo coth(-/3fiw) Im{x(w)}, 
7T 7 2 



(35) 



t P OO -1 

( p r) = ~ dojM 2 oj 2 coth{- f3hio)lm{x(uj)} 7 (36) 

7T J 2 



In this section the non-Markovian entanglement 
dynamics of a squeezed two mode state in the modi- 
fied common reservoir model is numerically analyzed. 
The initial covariance matrix is given by pip with 
elements (q 2 ) = (p 2 ) = icosh(2£) for j = 1,2 and 
({91,92}) = -{{pi,Pi}) = sinh(2£) (with m = 1, 
u>o = 1 and h — 1). The corresponding initial values 
of the variances of the normal coordinates therefore 
read (R 2 ) = ±e 2 «, (P 2 ) = e^, (x 2 ) = e^, and 
(Px) — i e2 ^ w hh minimal uncertainty (R 2 )(P^ / ) = 4 
and (a; 2 )(p 2 ) = \. 

Figures[T][4]give examples of the time evolution of a two 
mode squeezed states in a common reservoir, in form of 
the logarithmic negativity Ej^{t) ©(blue line), separa- 
bility function S(t) ([3]) (yellow line) and purity ji(t) (red 
line). While for small coupling parameter 7 and large 
cut-off frequency T the Markovian results would be re- 
produced, one can recognize deviations from the Marko- 
vian separability times T\ and T2 due to non-Markovian 
and strong coupling effects. The non-weak system- 
environment interaction accelerates decoherence with the 
consequence that the initial entanglement caused by the 
squeezing is lost faster (figflj. In contrast to the Born- 
Markovian results in ref. [1J| the initial entanglement is 
also lost if the critical value £ c is exceeded. Specific non- 
Markovian effects (which are relevant if the bath corre- 
lation time scale T _1 is comparable to the separability 
time scale) are shown in figlU Non-Markovian influences 
prolong the separability time r s and may lead to a Gaus- 
sian decay of the logarithmic negativity instead of an 
exponential decay. Furthermore, figH] shows that par- 
tial revivals of quantum correlations can occur after an 
initial loss of entanglement. This can be seen from the 
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FIG. 1: Time evolution of the logarithmic negativity Ej^(t) 
@ (blue line), separability function S(t) © (yellow line) and 
purity fi(t) (red line). The initial entanglement is lost when 
Em = and fi(t) is in the entanglement region ([9]) (shaded 
grey) respectively. Em = for t > t s where t s < Ti, T2 (here 
T\ ~ ujq 1 and T2 ~ 2.2WQ 1 ) due to the effect of non-weak 
interaction strength 7 = 0.2ljo- Parameters are |£| = 1 and 
F = 10lj , T = 3.5hw /k. 
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FIG. 2: Non-Markovian effects (reducing cut-off frequency to 
V > wo) prolong the separability time (compared to T > 
wo ) and may lead to a Gaussian decay of the logarithmic 
negativity. Partial revivals of the logarithmic negativity can 
occur. Here V — wo- Parameters £, 7, T and colors of the 
functions Ej^(t), S(t), fi(t) are chosen as in fig. 1. 



oscillations of the logarithmic negativity between positive 
values (entanglement) and zero (separability). The same 
conclusions can be drawn from the criterion for marginal 
and global purities. The two mode state is separable if 
the total purity /i(t) is in the region determined by eq. 
© (gray shaded area in figOH]) . 

The stationary state in figO] and [5] is separable in each 
case. However it is possible that quantum correlations 
exist in the stationary state. This becomes obvious from 
fig. [3] and 21 which show examples for low temperatures. 
Despite a loss of entanglement at short times t < t s quan- 
tum correlations between the two modes emerge again at 
later times and persist even in the steady state. Thus, 
in this model there is not just a separability or disentan- 
gling time r s but also a re-entangling time r e after which 
the two mode state remains entangled asymptotically. 
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Wot 

FIG. 3: At low temperature and strong coupling, the two- 
mode state can get re-entangled after a certain time r e with 
asymptotically constant quantum correlations. Parameters 
here are 7 = 1.5wo, T — IOljo and T = 10~ 3 hwo/k (with 
Markovian separability times T\,Ti — > 00). Symbols and col- 
ors are chosen as in Eg. 1. 
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FIG. 4: Increasing the coupling strength 7 reduces the separa- 
bility time T a as well as the re-entangling time r e (compared to 
fig. 3). Initial entanglement is destroyed faster while asymp- 
totic quantum correlations are generated sooner. Parameters 
here 7 = 2wo, F = lOwo and T = W~ 3 hwo/k. Symbols and 
colors are chosen as in Eg. 1. 



However it is difficult to find this characteristic time 
scale numerically due to the possible oscillations with 
partial revivals of the logarithmic negativity. Increasing 
the coupling strength 7 reduces the separability time r s 
as well as the re-entangling time r e as can be seen from 
a comparison between fig(3] and figdl The initial entan- 
glement resulting from the squeezing is destroyed faster 
whereas environment-induced asymptotic quantum cor- 
relations are generated sooner. The re-entanglement time 
scale is strongly affected by the system bath interaction 
strength 7 but the initial squeezing with parameter £ 
has little effect on the stationary quantum correlations, 
in contrast to the results in previous works [14L l3l| . Non- 
Markovian effects (induced by small cut-off frequency T) 
just extend the separability time scale r s but usually do 
not influence the re-entangling time r e . Entanglement 
of the two modes can exist in the stationary state if the 
bath temperature is below a critical value, that depends 



7 





0.25 





















. 2 


T(n\ 




\\Pl,P2}) 




. 15 








> 


0.1 










. 05 


Er x \ 


S 








. 05 





0.2 0.4 0.6 0.8 1 



kT/hu>o 

FIG. 5: Separability criterion S(oo) of eq. ([3]), logarithmic 
negativity Em (O and Er x — mmjj. — {ii 2 ){p^),0} com- 
pared to the mutual information I(p) (|40|) and the stationary 
correlations {{p 1 ,p2}), {{91,92}) for the two mode state in 
dependence of the temperature of the common reservoir. Pa- 
rameters 7 = 2wo and V = 10wo- 




FIG. 6: Stationary entanglement measured by the function 
Er x = max{^ — (R 2 ){p x ) , 0} fx 10~ 4 j of a two-mode state 
in a common reservoir in dependence of coupling strength 7 
and cut-off frequency F at temperature T = 0.25hu>o/k. 



on the coupling parameter 7 and the cut-off frequency T. 
In the following we derive a simple separability criterion 
for the stationary state. Since the normal coordinates 
R = (qi + ©)/2 and p^ = (pi - p 2 )/2 are EPR-like 
operators we can use product version |(7J) of the EPR- 
separability criterion in the form 

(R 2 ){pl)>^ O ((<Zi+92) 2 )((pi-P2) 2 )>1 (38) 

in order to derive a necessary condition for the existence 
of quantum correlations in the stationary state. Since in 
this case ({91,52}) < 0, {{pi,p 2 }) > and the stationary 
covariance matrix of the modified QBM model is given 
in symmetric standard form, the pro duct criterion ([7]) 
is equivalent to the PPT-criterion [51| and thus detects 
the same set of entangled covariance matrices. Using the 
stationary correlations in eq. and from (|57| this 

separability criterion can be rewritten in the form 

ui f°° 1 1 

— / doj{2n+l)coth(-Phaj)lm{x(uj)} > (39) 

7T JO 2 4 

with mean occupation number n = {e^ hu) ° — of the 
relative coordinate. Figure [5] compares this criterion 
in form of the function Er x = max{^ — (i? 2 )(p^),0} 
with the partial transposition separability criterion ([3]) 
and the logarithmic negativity (0 in dependence of the 
bath temperature. All these criteria give the same critical 
temperature. Also plotted are the stationary correlations 
({<7i, 92}) and ({pi,P2}) as well as the mutual information 
of a two mode Gaussian state which is defined by 



I(p) = S„(pi) + S v (p 2 ) - S v (p), 



(40) 



where S v (p) is the von Neumann entropy of the total 
two mode state and S v (pi), S v (p2) are the von Neumann 
entropies of the reduced states p\ — Tr 2 [p] and p 2 = 
Tri [p] . The von Neumann entropy of the reduced systems 



is equivalent to that of a single mode state and is given 
by 



S v (pj) 



-In 



2fij 



1 ~ N In 1 + Mj 
2pj 1 — Hj 1 + p 3 



(41) 



while the total entropy can be derived from the symplec- 
tic eigenvalues v± of V and reads [13, HH : 



S v (p) = fiy-) + f(y+) 



(42) 



with f{v±) - {v± + \) hx{v ± + i) - { V± - i) \n{v± - i). 
The mutual information is a measure of the amount of 
quantum and classical correlations [54[ and decays more 
slowly with temperature than the logarithmic negativity. 
Above a critical temperature (where the logarithmic neg- 
ativity becomes zero) only classical correlations remain. 
The critical temperature depends on the spectral density 
of the bath and on the response function of the subsys- 
tem. The existence of stationary quantum correlations in 
dependence of the coupling constant 7 and the cut-off fre- 
quency r is illustrated by figure[6] If the system-reservoir 
interaction constant becomes strong enough at a given 
temperature quantum correlations are still present in the 
stationary two-mode-state. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the dynamics of two 
mode squeezed states in a common thermal reservoir 
within an extended quantum Brownian motion model. 
Contrary to the two-reservoir model a common bath can 
provide an indirect coupling between the two subsystems 
and therefore a mechanism to correlate them. Thus, 
the coupling to the environment induces not only 
decoherence leading to separability but also can generate 
correlations resulting in asymptotic entanglement at 
low temperatures. The influence of non-Markovian and 
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strong coupling effects on separability times scales was 
compared to previous results for Markovian dynamics. 
Non-Markovian effects increase the separability time 
but have little effect on the re-entangling time. The 
re-entangling time depends first of all on the system bath 
interaction strength. Increasing this coupling constant 
enforces entanglement generation and therefore reduces 
the re-entangling time. However, it also enhances 
decoherence and thus weakens the initial quantum 
correlations reducing the separability time. For given 
parameters of the bath spectral density asymptotic en- 
tanglement is only present below a critical temperature. 
From the separability criterion for EPR-like operators a 
separability criterion was derived, which depends on the 



bath temperature and the response function of the open 
quantum system. Summarizing, the extended quantum 
Brownian motion model of a two-mode continuous vari- 
able system in a common reservoir provides an example 
of the case where environment-induced entanglement 
generation counteracts environment-induced decoher- 
ence. In view of the theory of quantum information 
with continuous variable systems these results may 
be relevant for experimental applications of correlated 
quantum channels. 

Acknowledgement - We would like to thank J. Eisert 
for valuable comments on the subject of this paper. 



S. L. Braunstein and A. K. Pati, Quantum Information [22 
Theory with Continuous Variables (Kluwer, Dordrecht, 
2003). [23 
S. Braunstein and P. van Look, Rev. Mod. Phys. 77, 513 
(2005). [24 
N. J. Cerf, G. Leuchs, and E. S. Polzik, Quantum Infor- 
mation with Continuous Variables of Atoms and Light [25 
(Imperial College Press, London, 2007). 

A. Furusawa, J. L. Sorensen, S. L. Braunstein, C. A. [26 
Fuchs, H. J. Kimble, and E. S. Polzik, Science 282, 706 [27 
(1998). 

S. L. Braunstein and H. J. Kimble, Phys. Rev. Lett. 80, [28 
869 (1998). 

R. Simon, Phys. Rev. Lett. 84, 2726 (2000). [29 
L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 84, 2722 (2000). [30 
G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 [31 
(2002). [32 

G. Giedke, M. M. Wolf, O. Krger, R. F. Werner, and J. I. 
Cirac, Phys. Rev. Lett. 91, 107901 (2003). [33 
Z. Y. Ou, S. F. Pereira, H. J. Kimble, and K. C. Peng, 
Phys. Rev. Lett. 68, 3663 (1992). [34 
T. C. Zhang, K. W. Goh, C. W. Chou, P. Lodahl, and 

H. J. Kimble, Phys. Rev. A 67, 033802 (2003). [35 
H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, Oxford, [36 
2003). [37 
T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, [38 
G. Schon, and W. Zwerger, Quantum Transport and Dis- [39 
sipation (Wiley- VCH, Weinheim, 1998). [40 
J. S. Prauzner-Bechcicki, J. Phys. A 37, L173 (2004). 

D. Braun, Phys. Rev. Lett. 89, 277901 (2002). [41 
A. Serafini, F. Illuminati, G. A. M. Paris, and 
S. De Siena, Phys. Rev. A 69, 022318 (2004). [42 
S. Olivares, M. G. A. Paris, and A. R. Rossi, Phys. Lett. 
A 319, 32 (2003). [43 
J. J. Halliwell, J. Perez-Mercader, and W. H. Zurek, 
Physical Origins of Time Asymmetry (Cambrige Univ. [44 
Press, Cambrige, 1994). 
M. S. Kim, J. Lee, D. Ahn, and P. L. Knight, Phys. Rev. [45 
A 65, 040101 (2002). [46 
X. X. Yi, C. S. Yu, L. Zhou, and H. S. Song, Phys. Rev. [47 
A 68, 052304 (2003). 

F. Benatti and R. Floreanini, J. Phys. A 39, 2689 (2006). [48] 



F. Benatti and R. Floreanini, Int. J. Quant. Inf. 4, 395 
(2006). 

M. B. Plenio and S. F. Huelga, Phys. Rev. Lett. 88, 
197901 (2002). 

A. K. Rajagopal and R. W. Rendell, Phys. Rev. A 63, 
022116 (2001). 

D. F. Walls and G. J. Milburn, Quantum Optics 

(Springer, Berlin Heidelberg, 1994). 

M. Ban, J. Phys. A: Math. Gen. 39, 1927 (2006). 

H. McAneney, J. Lee, D. Ahn, and M. S. Kim, J. Mod. 

Opt. 52, 935 (2005). 

S. Maniscalco, S. Olivares, and M. G. A. Paris, Phys. 
Rev. A 75, 062119 (2007). 

J. H. An, M. Feng, and W. M. Zhang, quant- 
ph/0705.2472 (2007). 

K.-L. Liu and H.-S. Goan, quant-ph/0706.0996 (2007). 
J. H. An and W. M. Zhang, quant-ph/0707.2278 (2007). 

A. O. Caldeira and A. J. Leggett, Physica A 121, 587 
(1983). 

G. W. Ford, J. T. Lewis, and R. F. O'Connell, J. Stat. 
Phys. 53, 439 (1988). 

B. L. Hu, J. P. Paz, and Y. Zhang, Phys. Rev. D 45, 
2843 (1992). 

G. Adesso, S. A., and F. Illuminati, Phys. Rev. A 70, 
022318 (2004). 

S. M. Tan, Phys. Rev. A 60, 2752 (1999). 

M. D. Reid, Phys. Rev. A 40, 913 (1989). 

V. Giovannetti, Phys. Rev. A 67, 022320 (2003). 

P. Hyllus and J. Eisert, New J. Phys. 8, 51 (2006). 

G. Adesso, A. Serafini, and F. Illuminati, Phys. Rev 

Lett. 92, 087901 (2004). 

G. Adesso, A. Serafini, and F. Illuminati, Open Syst. Inf. 
Dyn. 12, 189 (2005). 

A. O. Caldeira and A. J. Leggett, Ann. Phys. (NY.) 149 
374 (1983). 

A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46 
211 (1981). 

G. W. Ford, M. Kac, and P. Mazur, J. Math. Phys. 6 
504 (1965). 

R. Karrlein and H. Grabert, Phys. Rev. E 55, 153 (1997) 
F. Haake and R. Reibold, Phys. Rev. A 32, 2462 (1985) 
J. E. und M.B. Plenio, Phys. Rev. Lett. 89, 137902 
(2002). 

C. -H. Chou, T. Yu, and B. L. Hu, quant-ph/0703088 



9 



(2007). 

[49] C. Horhammer and H. Biittner, J. Phys. A: Math. Gen. 

38, 7325 (2005). 
[50] O. S. Duarte and A. O. Caldoira, Phys. Rev. Lett. 97, 

250601 (2006). 

[51] M. S. Kim and J. Lee, Phys. Rev. A 66, 030301 (2002). 
[52] A. Serafini, F. Illuminati, and S. De Siena, J. Phys. B 



37, L21 (2004). 
[53] A. S. Holevo and R. F. Werner, Phys. Rev. A 63, 3032312 

(2001) . 

[54] H. Ollivier and W. C. Zurek, Phys. Rev. Lett. 88, 017901 

(2002) . 



